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OSCULATORY INTERPOLATION FORMULAS. 
By C. H. Fokstth, Ann Arbor, Mich. 



The best formulas of interpolation, making use of finite 
differences, are Newton's, Everett's, Stirling's, and Bessel's 
which are as follows: 

Newton's formula: 

, . , x(x — 1).„ . x(x — l)(x — 2) ., . 
u x = M + xAu H — y A 2 M -1 — gy A M o +• • • • 

Everett's formula: 
u x = f(g)u +f(x)ui 

, (z + l) (3 > , (z + 2) (5 > , 
where f(x)=x+ v ^ y 5 2 + v -^ y S 4 

/ . (s+l)s(a:-l) , 2 , (x+2)(x+l)x(x-l)(x-2) d 

= { X+ 3! S + 5! 5 

and £ is the distance of u x measured behind «i. 

Stirling's formula: 

, n 2 , n(n 2 — 1)-, , n 2 (n 2 — 1) .. 
u x =Uo+nbu a + — S 2 Uo-\ r^ 5 3 W(H jj ^ 3 M ° 

+ B ' (n, - 1 K ) , (B '-' 4) ^+.- 

5! 



Bessel's formula: 

n(nj-i)- n(n-l)(n -f), 
2! * i+ " 3! 



, /mm. — j. / -„ . ivy rv — x i y iv — ^ / . 

u x = u + nbu>+ ol d*u t + — 57 -&u h 



n(n 2 -l)(n-2) ^ | n(n 2 - l)(n-2)(n-f) ^ + 



4! 5 5! 

Thus Newton's formula is based upon ordinary differences 
and the three others upon central differences whose notation 
is explained by the following table: 
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«-» 

8u- t 

Uo 5 2 m S 4 u 

5mj 8 s u± tf"u§ 

«! 5 2 wi 5 4 wi 

8u% 5 3 Mg 

MS 

We shall use the central difference symbols but with a dash 
over them to designate the arithmetic averages of the differ- 
ences of the same order on both sides of them. Thus 
S«o = §(8wj+S»-j), and so on for the Best. We believe that this 
latter idea is an improvement over that of Mr. Sheppard,* 
who contrived the scheme of representing central differences 
shown in the above table, in which a separate character n is 
introduced instead of the dash which we suggest. We con- 
sider the use of the character /i in such a way to be very 
misleading. 

The purpose of this paper is threefold. First we wish to give 
the results of our derivations of the osculatory formulas based 
upon Stirling's and Bessel's formulas corresponding to those 
which have been derived based upon Newton's t and Everett's J 
formulas, all to be used for interpolating four values between 
two of a series of equidistant values by fifth differences. Sec- 
ond, we wish to bring together all four osculatory formulas 
for purposes of comparison and future reference. Third, we 
wish to offer for future reference formulas for computing fifth 
leading differences to be used in connection with the applica- 
tions of the four osculatory interpolation formulas, to simplify 
the computation when several successive intervals are to be 
considered. 

As is well known, all interpolations within a series of suc- 
cessive intervals are performed through the use of a series of 

♦Proceedings of the London Mathematical Society, Vol. XXXI, p. 149. 
t Journal of the Institute of Actuaries, Vol. 22., p. 270. 
JJ. I. A., Vol. 42, p. 369. 
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what might be called partial interpolation curves. In os- 
cillatory interpolation, the "joints" of these separate inter- 
polation curves are welded together in a more satisfactory 
manner by requiring that not only the slopes but also the 
curvatures of any two interpolation curves to be the same 
at their point or points of intersection. Thus the two inter- 
polation curves, one through the points u-2, u-i, Wo, Ui, and Uz, 
the other through w-i, w , Mi, u 2 , and u% are required to have 
equal slopes and equal curvatures at the two points. m and 
wi for purposes of interpolating any number of values (four, 
in this paper) between w and u\. This modification tends to 
remove the discontinuities usually found at the intersections 
of two distinct interpolation curves. 

This method of interpolation has been termed "osculatory" 
interpolation because the successive pairs of partial interpola- 
tion curves are thus seen to have a common osculating circle 
at their points of intersection. It was first devised by Sprague* 
and later developed and applied by Karup, t King, J Buchanan,! 
and others. || The method or rather modification has already 
been applied to Newton's* and Everett's§ formulas, the latter 
of which involves only even central differences of each of the 
middle pairs of the series of values between which the inter- 
polations are to be made. The central difference formulas, 
Stirling's f and Bessel's, are so widely used that it seemed well 
worth while to obtain the osculatory formulas for them as 
well as for Newton's and Everett's formulas. 

The process of deriving the osculatory formulas wherein 
analytically first derivatives are considered in lieu of slopes 
and second derivatives in lieu of curvatures has been explained 
so well in connection with the derivation of Sprague's formula 
(based upon Newton's formula) and Buchanan's formula 
(based upon Everett's formula), and especially in an independ- 
ent note by Lidstone**, that we shall give simply the final 

* Journal of the Institute of Actuaries, Vol. 22, p. 270. 
t Trans. Second International Actuarial Congress, p. 78. 
JJ. I. A., Vol. 42, p. 225. 

I J. I. A., VoL 42, p. 369. 

|| Quarterly Pub. Amer. Stat. Assoc, June, 1910. 
Tie Record. Amer. Inst, of Actuaries, June, 1911. 

II Text-book. Institute of Actuaries, p. 447. 
** J. Inst. Act., Vol. 42, p. 394. 
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equations of the two curves. The one based upon Stirling's 
formula becomes 

u x = Uo+ X5u + ^■5 2 Uo+ x —^-. — -5 3 Wo+- — —. — -8% 

3. 4. (i) 



4. (— _ — _L ^\ i^f? 

+ \20 4 + 5/ 3! 



The equation based upon BesseFs formula becomes 
u x =u + x8u i + ^^lhi i + (tf-%x* + ^ S ^ (2) 



(x* - 2a; 3 - x 2 + 2»)-. , (-a; 6 + 2 a; 4 -a; 3 ) 
" — 4| -«S + - 4] " 



If we substitute successively x = %, i, ■■■, | in equations 
(1) and (2), and difference each of the sets of six results five 
times we obtain the two desired sets of leading differences 
for interpolating four values between «o and «i. 
The osculatory leading differences based upon Stirling's 
formula are: 



(A) 



The osculatory leading differences based upon BesseFs 
formula are: 



(1) 


5Mo,! S 2 u . S s u S 4 u a 5% 
5 2 5 2 5 3 5 4 5 4 
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" 1 1 " 
T" 2 
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" + " 
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(B) 



The osculatory leading differences derived by Buchanan, 
based upon Everett's -formula are: 
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(1) 


Mo „ 5 2 M . „ 


(2) 


4 " -7 


(3) 


-1 " +2 


(4) 


3 


(5) 


-5 



a^o mi . a 2 Mi a 4 Mi 

5 4 + 5 5 3 + 5 4 

+ 1 " +3 " 

+1 " -3 " 

-2 " 

+ 5 " 



(C) 



The osculatory leading differences derived by Sprague, 
based upon Newton's formula are: 



,.,, AMo A 2 m , ,, A 3 m n A 4 M . A S M 

(1) ^ _ + 8_ 5^ +11 "5^~ 11 ^" + "s 4 " 

(2) 1 " +6 " +1 " +3 " 

(3) 1 " +4 " -3 " 

(4) 1 " -2 " 

(5) 5 " 



(D) 



To apply these different sets of leading differences, the 
observed series of values are differenced as usual, the differ- 
ences averaged wherever necessary in connection with the 
central differences and finally each difference divided by the 
appropriate power of 5 (or multiplied by the corresponding 
decimal) . The leading differences themselves are then formed 
according to the above tables. Each group of six original 
values will have its own set of differences for interpolating 
four values between the inmost pair, but if a series of succes- 
sive intervals are to be interpolated the best plan is to inter- 
polate the first interval as outlined above and then compute 
merely the appropriate fifth leading differences for the suc- 
ceeding intervals such that by continued summation or addi- 
tion to the preceding leading differences of lower order (in 
each case) the required interpolations will be obtained. The 
possibility of this plan has been pointed out before* and we 
shall merely carry it out as Buchanan has done in connection 
with his treatment of Everett's formula. Not being able to 
find the formula for computing these fifth differences in connec- 
tion with Sprague's formula, we have derived these formulas 
for this case as well as for the osculatory formulas based upon 

*J. I. A. Vol. 42, p. 285. 
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Stirling's and Bessel's and derived in this paper. The method 
of derivation is outlined by Buchanan in his paper. These 
formulas or corrections of fifth differences for all four oscu- 
latory formulas are as follows: 





Newton's 


Everett's 


Stirling's 


Bessel's 
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The greatest disadvantage found in connection with the 
use of the above fifth differences is that usually the work 
should be carried to a greater number of decimals than if 
the above leading differences were used. If an appropriate 
computing machine is available this criticism is of no import- 
ance. As a general rule the work in connection with the 
leading differences should be carried to two more decimals 
than are to be finally retained and the work of the fifth 
differences to three or the maximum four. 

One great advantage in the use of the above osculatory 
formulas is the accompanying check upon the computation, 
for the sets of leading differences not only lead to interpola- 
tions of the four values between a pair of intermediate original 
values but will also reproduce these original values and since 
the process of continued summation causes any errors in the 
computation to accumulate, such errors are readily caught. 

Without going into the details, it should perhaps be pointed 
out that Bessel's and Everett's and hence the corresponding 
osculatory formulas lead to about the same degree of accuracy 
but are generally considered to be the most accurate of the 
four formulas considered in this paper. However, Everett's 
formula and the corresponding osculatory formula involve 
considerably less as well as simpler computation than Bessel's 
formula and its corresponding osculatory formula. The oscula- 
tory formula based upon Stirling's formula is very little less ac- 
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curate than those based upon Everett's and Bessel's formulas 
and is remarkably easy to apply although scarcely as easy as 
Buchanan's (that based upon Everett's). There is scarcely 
any justification in the further use of Sprague's osculatory 
formula since the derivation of the formulas for the other 
osculatory interpolations, leading as it does to results nowise 
as accurate and involving much computation. 

In conclusion the following example is given to show the 
use of the above sets of leading differences. The following 
values of u are arbitrarily chosen and the corresponding 
differences found in the usual way as shown. 

u-%= 4673423 

13462 
w-i =4686885 6321 

19783 384 

u =4706668 6705 65 

26488 449 13 

mi =4733156 7154 78 

33642 527 

m 2 =4766798 7681 

41323 
« 3 =4809121 

Here we have, for example, 5m = £(19783+26488) 
The results after applying the four sets of leading differences 
are: 

By (A) By (B) By (C) By (D) 

Mo =4706668 4706668 4706668 4706668 

4711406 4711397 4711397 4711416 

4716418 4716401 4716401 4716435 

4721710 4721691 4721691 4721729 

4727287 4727273 4727273 4727301 

Mi =4733156 4733156 4733156 4733156 



(C)-(A) (C)-(B) (C)-(D) 
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